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Abstract 

^ : 

We evaluate the effect of the bottom-quark mass on resummed transverse momentum (qx) distri- 
butions of supersymmetric Higgs bosons at the Tevatron and LHC. The mass of the bottom quark 
acts as a non-negligible momentum scale at small qr and affects resummation of soft and collinear 

•rH , 

^ . radiation in this region. The improved treatment of the fe-quark mass and kinematical effects leads 

h : 

to observable modifications in the resummed predictions for both colliders. 
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I. INTRODUCTION 



Understanding the nature of electroweak symmetry breaking is the central challenge 
for high-energy physics. The search for Higgs bosons, assumed to be responsible for the 
generation of gauge-boson and fermion masses, is the primary task for the existing and 
future colliders. 

The Higgs sector may be represented by one complex scalar doublet, as it is economically 
realized in the Standard Model (SM), or by two or more doublets, as it takes place in 
the Minimal Supersymmetric Standard Model (MSSM) and its extensions. An important 
feature of MSSM is that, for large values of tan/5, the Yukawa couplings of the 6-quarks to 
the neutral Higgs bosons TC (where TC — h, H, or A) are strongly enhanced compared to the 
SM bbHsM Yukawa coupling. Consequently, production of supersymmetric Higgs bosons in 
bb fusion can have a large cross section in supersymmetric extensions of the Standard Model 

D □ □ □ □ 

The partonic processes contributing to the inclusive Higgs boson production with en- 
hanced bbTC coupling are represented by (a) bb — > TC; (b) gb — > TCb; and (c) gg — > 
bbTC scattering. The three processes (a,b,c) all give rise to the same hadronic final 
states, with two -B-mesons appearing in different, but overlapping, regions of phase 
space. The distinction between the three processes depends very much on the factoriza- 
tion scheme adopted for the QCD calculation, as has been recently reviewed in Ref. 
The TCbb processes have been extensively studied recently in SM and MSSM scenar- 
ios Q, S fl H fl 0, H Ji4 , H Q, H Q ■ 

As shown in Refs. |l9l. l2fl|. the correct model for the transverse momentum distribution 
of the Higgs boson is crucial for unambiguous reconstruction of the Higgs boson mass in 
the TC — > tt decay channel. It is also important for discriminating the signal events from 
the backgrounds by examining the qx distribution of the Higgs boson in TCbb associated 
production, followed by TC — ► bb decay [2j]]. The transverse momentum (gy) distributions of 
Higgs bosons may be sensitive to the mass of the bottom quark when qx is comparable 
to m&. In Refs. [3, 13| , we study the effect of the initial-state multiple parton radiation 
and heavy-quark masses on the transverse momentum distribution in the bb — > TC process. 
Here we summarize the results of those two papers. 
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II. TRANSVERSE MOMENTUM RESUMMATION FOR MASSIVE QUARKS 



The resummed differential cross section for inclusive production of Higgs bosons in scat- 
tering of initial-state hadrons A and B takes the form (2J| 

da r°° bdb 



dQ 2 dydq^ 



f 00 bdb 

/ — — </o(<?rb) W(b, Q,x A ,x B ,m b ) + Y(q T ,Q,y,m b ) } (1) 
Jo 27T 



where y is the rapidity of the Higgs boson, xa,b = Qe ±y /vS are the Born-level partonic 
momentum fractions, S is the square of the center-of-mass energy of the collider, and Jo(gTb) 
is the Bessel function. The resummed form factor W is given in impact parameter (b) space 
and factorizes as 

W(b, Q, x A , x B , m b ) = ^J2 a fk e~ s{h ' Q ' mb) V j/A (x A: b, m b ) V k / B (x B , b, m 6 ), (2) 

where the summation is performed over the relevant parton flavors j and k. Here, is a 
product of the Born-level prefactors, e - s ( h <Q< m b) [ s an exponential of the Sudakov integral 

[ q2 dp, 2 r _ (Q 2 \ _ ] 

S(b,Q,m b ) = / —^A(a s ([i),m b )\n( — )+B(a s ([i),m b ) , (3) 

with b = 2e~ lE ~ 1.123, and Vj/ A (x,b,rn b ) are the b-dependent parton distributions for 
finding a parton of type j in the hadron A. In the perturbative region (b 2 A^), the 
distributions Vj/ A (x,b,m b ) factorize as 

P i/j4 (x,b,m 6 )| b2<<A _ 2 = ^2 -rCj/a(x/£,b,m b ,fi F )f a/A (£,fi F ) (4) 

QCD j Jx S 

into a convolutions of the Wilson coefficient functions Cj/ a (x, b, m b , hf) and fcr-integrated 
parton distributions f a /A(£, A*f)- The Sudakov exponential and b-dependent parton den- 
sities resum contributions from soft and collinear multi-parton radiation, respectively. 
Y = PERT — ASY is the difference between the finite-order cross section (PERT) and 
its asymptotic expansion in the small-g^ limit (ASY). 

The Higgs cross sections depend on the mass m b of the bottom quark. The distributions 
Vj/ A (x,b,m b ) for the heavy quarks (j = c,b) cannot be reliably evaluated at all impact 
parameters if a conventional factorization scheme, such as the zero-mass variable-flavor 
number (ZM-VFN, or massless) scheme, is used. The reason is that m b acts as an additional 
large momentum scale, which, depending on the value of b, introduces large logarithms 
\n n (m b b) or non-negligible terms cx {m b b) n . The situation encountered here is reminiscent 




Figure 1: The bottom-quark distributions V b / p (x,b,mb) in the proton vs. the impact parameter b. 
The solid and dashed curves correspond to the S-ACOT and massless ("ZM-VFN") factorization 
schemes, respectively. 

of the heavy-quark contributions to the DIS structure functions Fi(x,Q 2 ), which are not 



adequately described by the conventiona 
mentum transfers Q 2 (see, for instance, 



factorization schemes at either small or large mo- 
To work around this complication, Ref. [3| 
proposed to formulate the CSS formalism in a general-mass variable flavor number (GM- 
VFN) scheme j^, which correctly evaluates the heavy-quark mass effects at all momentum 
scales. Among all GM-VFN factorization schemes, the S-ACOT scheme [3, 13| was found 
to be well-suited for the efficient calculation of the CSS resummed cross sections. In partic- 
ular, in this heavy-quark CSS (CSS-HQ) formalism Q| the dependence on m& is dropped 
in all 0(a s ) terms in Eq. (JTJ) except for Vb/A{%, b, m&). 

The dependence of the bottom-quark parton density Vb/ P (x, b, m&) on the impact param- 
eter is shown in Fig. [TJ The ZM-VFN parton density Vb/ P (x, b, m&) is not properly defined 
below the threshold \ip = nib (or above b = bo /nib). It was continued to large b in the pre- 
vious calculations using an effective "ZM-VFN" approximation described in Ref. The 
S-ACOT parton density Vb/ P (x,b,nib) is well-defined at all b. It reduces to the ZM-VFN 
result at b <C b Q /m b and is strongly suppressed at b 3> b /m b . The suppression is caused 
by the decoupling of the heavy quarks in the parton densities at /j,p much smaller than nib 
(b much larger than bo /nib). Consequently the impact of the non-perturbative contributions 
from b > 1 GeV -1 is reduced in the heavy-quark channels compared to the light-quark 
channels. 



4 



pp (bb) ^HX (VS = 1960 GeV) 



pp (bb) ->HX (Vs = 14 TeV) 



CTEQ5HQ 1 
M„ = 120 GeV 




- 


CTEQ5HQ1 : 




IVljj — 1 L\) uc V 


\ j Massive (S-ACOT) 




" / Massless ("ZM-VFN") 





15 20 25 30 
q T [GeV] 



45 



10 



15 20 25 
q T [GeV] 



30 35 40 45 



(a) 



(b) 



Figure 2: Transverse momentum distribution of on-shell Higgs bosons in the bb —* TL channel 
at (a) the Tevatron and (b) LHC. The solid (red) lines show the qx distribution in the massive 
(S-ACOT) scheme. The dashed (black) lines show the distribution in the massless ("ZM-VFN") 
scheme. The numerical calculation was performed using the programs Legacy and ResBos Hi 
with the CTEQ5HQ1 parton distribution functions [31]. The bottom quark mass is taken to be 
m b = 4.5 GeV. 



The massless ("ZM-VFN") calculation therefore underestimates the true behavior at 
b > 0.1 GeV -1 and small qr- This effect can be seen in Fig. which displays dojdqx for 
bb — > ri boson production at (a) the Tevatron and (b) LHC. 1 At the Tevatron, the qr max- 
imum shifts in the "ZM-VFN" approximation to larger qr by about 2 GeV out of 11.7 GeV 
(about 17%). For a Higgs mass M H = 200 GeV, the maximum of da/dqx shifts by about 
1.9 GeV out of 12.7 GeV. At the LHC, the difference between the "ZM-VFN" and S-ACOT 
calculations is smaller compared to the Tevatron, because the influence of the b > 0.1 GeV -1 
region is reduced at smaller momentum fractions x probed at the LHC [^2|. The maximum 
of the qr distribution shifts in the "ZM-VFN" approximation by about 1.3 GeV (9% out of 
14.1 GeV) to larger qx- The results for other Higgs masses can be found in Ref. 



1 Fig. does not specify the overall normalization of qr distributions. It is valid for both Standard Model 
and supersymmetric Higgs bosons, since at leading order the supersymmetric result can be obtained by 
rescaling the Standard Model bbHsM coupling: 9^1^ ^ A\ = {~ sina, cosa, sin/?75}gj^|/cos/3. The net 
effect of mb on qx distributions will be the same for the SM and MSSM Higgs bosons, up to an overall 
normalization constant. 
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Figure 3: qr distributions for production of 100 GeV CP-odd Higgs bosons A via bb fusion in the 
Tevatron Run-2. The solid and dashed curves correspond to the lowest-order W-term W(l, 1, 0) 
(with functions A(a s (p)) and B(a s (p,)) evaluated at 0(a s )) and PYTHIA. 

III. NUMERICAL COMPARISON WITH PYTHIA 

The full qr dependence of the bb — > H process is affected by constraints on phase space 
available for QCD radiation (less relevant at small qr). We illustrate the interplay of various 
effects by comparing the CSS-HQ resummation to the PYTHIA Monte Carlo program 
We focus on production of the CP-odd Higgs particle A for tan (3 — 50 (predictions for the 
other Higgs bosons can be obtained by rescaling the bbA coupling). 

As compared to the CSS-HQ formalism, the PYTHIA calculation does not include con- 
tributions generated from the C-functions and Y-term, and it evaluates the soft parton 
contributions at 0(a s ). Therefore, we start by comparing the PYTHIA qr distribution to 
the resummed W-tevm W(l, 1, 0) in Eq. (IT}, with the functions A, B, and C in Eqs. Q, © 
being evaluated at orders a s , a s , and a°, respectively. The orders of a s in A, £>, and C are 
shown as the arguments of W(l, 1, 0). 

It is evident from Fig. El that the shapes of W(1,1,0) and PYTHIA qr distribution 
are very different, though the integrated rates (i.e., the areas under the two curves) are 
about the same. The qr distribution from PYTHIA is narrower and peaks at lower qr than 



6 



W(l, 1, 0). The large discrepancy between the two curves is in contrast to the case of W and 
Z production via light-quark scattering, where the above two calculations predict similar, 
though not identical, qr distributions [29]. 

A closer examination reveals that additional features must be implemented in the re- 
summed cross section in order to reliably describe the qr distributions of Higgs bosons 
produced via bb fusion. 

• The kinematical effects account for a large part of the disparity between W(l,l,0) 
and PYTHIA. The bottom-quark PDF is a rapidly decreasing function of x in the 
probed range of x. Consequently, approximations for the true partonic kinematics 
(especially those made for the light-cone momentum fractions x) may have a strong 
impact on the rate of bb scattering. This feature should be contrasted to the behavior 
of the light-quark PDF's in W and Z production, which include a substantial valence 
component and vary slower with x. As a result, the kinematical approximations are 
less consequential in the W and Z case. 

When PYTHIA generates QCD radiation, the kinematical distributions of the 
final-state particles, including the quarks and gluons from the QCD showering, are 
modified to satisfy energy-momentum conservation at each stage of the showering. 
In the resummation calculation, information about the exact parton kinematics is 
included in the finite-order term (PERT). The resummed cross section is therefore 
expected to be closer to PYTHIA once the 0(a s ) finite term, PERT(1)-ASY(1), is 
included. In the W(l, 1,0) calculation, the emitted gluons are assumed not to carry 
any momentum at all in the soft limit. To compensate for small, but nonzero energy 
of the soft gluon emissions, we introduce a "kinematical correction" (KC) in the W 
and ASY terms. This correction modifies the minimal values of partonic momentum 
fractions xa and xb in order to account for reduction of phase space available for 
collinear QCD radiation at large qr- 

• The lowest-order cross section W(l, 1, 0) does not evaluate effects of the bottom-quark 
mass, which is first included in the C-function of order a s . Also, additional, though 
not complete, 0(a 2 s ) contributions arise in the Sudakov form factors inside PYTHIA 
when the next-to-leading order PDF's are used. To account for both features, we 
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(a) (b) 

Figure 4: Comparison of qr distributions predicted by TOT(l), PERT(l) and PYTHIA, for Higgs 
boson produced via bb fusion at (a) the Tevatron Run-2 and (b) LHC, for Ma =100 and 300 GeV 
respectively. 

evaluate the W term at one higher order (2,2,1) and include the dependence using 
the CSS-HQ scheme. 

Thus, our full prediction TOT(l) is obtained by adding ^ (2,2,1) (evaluated 

in the CSS-HQ formalism with the kinematical correction) and PERT(l), and subtracting 
ASY KC (1). It is shown for M A = 100 GeV at the Tevatron in Fig. H^a) and M A = 300 GeV 
at the LHC in Fig. Efb). TOT(l) (solid line) is compared to the fixed-order prediction 
PERT(l) (dashed) and the PYTHIA prediction (dot-dashed). As one can see, the results 
for Tevatron and LHC are qualitatively similar. TOT(l) is closer to the PYTHIA prediction 
than W(l, 1, 0), though the two distributions are not identical. The PYTHIA qx distribution 
peaks at lower qx than TOT(l). In the large qr region, the TOT(l) rate is larger than the 
PYTHIA rate. 

Finally, Fig. El shows the integrated cross section as a function of the minimal qx in the 
calculation for the Tevatron (left) and LHC (right). This is another way to illustrate the 
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Figure 5: Comparison of the integrated rates, deduced form Fig.QJ as a function of the minimal qr 
value taken in the integration over qr at the Tevatron Run-2 (left) and LHC (right) for Ma =100 
and 300 GeV, respectively. 

differences in the shapes of qr distributions obtained in the resummation, fixed-order, and 
PYTHIA calculations. 



IV. CONCLUSION 

Multiple parton radiation in 6-quark scattering is conspicuously sensitive to effects of 
large bottom-quark mass m& and phase-space constraints on collinear emissions. Both m\, 
dependence and phase-space dependence tangibly modify the shape of Higgs qr distributions 



in the bb — > 7i processes. The two types of effects were consistent 
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y implemented within the 
291 , realized in a massive 



CSS resummation formalism for heavy-quark scattering 
(GM-VFN) factorization scheme. These corrections act on different qr regions. When the 
dependence on ra& is taken into account, the position of the peak in the da/dqr distribution 
shifts to a lower qr value, leaving the rate at large qr essentially unchanged. The kinematical 
correction is effective in the high-g^ region, where it largely reduces the Higgs production 
rate. As a result, we obtain an improved prediction for the full qr spectrum of Higgs bosons, 
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an important piece of information needed for the future Higgs searches. 
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